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The B — > K* p, K*uj decays are useful to determine the CKM angle 03 = 7. Their polarization 
fractions are also interesting since the polarization puzzle of the B — ► <f)K* decay. We study these 
decays in the perturbative QCD approach based on kx factorization. After calculating of the 
non-factorizable and annihilation type contributions, in addition to the conventional factorizable 
contributions, we find that the contributions from the annihilation diagrams are crucial. They 
give dominant contribution to the strong phases and suppress the longitudinal polarizations. Our 
results agree with the current existing data. We also predict a sizable direct CP asymmetries in 
B + — ► K* + p°, B° — ► K* + p~, and B + — > K* + u> decays, which can be tested by the oncoming 
measurements in the B factory experiments. 
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Table 1: Longitudinal Polarization fractions of some B — > VV moc 


es. 


Process 


Belle 


Babar 


QCDF [14, 24] 


QCDF±FSI[17] 


5 U -► (j)K*° 


0.45 ± 0.05 ± 0.02 


0.52 ±0.05 ±0.02 


0.91 


43+ u - i3 


B + -> 


0.52 ±0.8 ±0.03 


0.46 ±0.12 ±0.03 


0.91 


43+ u - i3 


5+ -► p u A^*+ 




°- 96 -ai5 ± °- 04 


0.94 


4Q+ U - n 


5+ -► p+£T* 


0.43 ±0.11+^ 


0.79 ± 0.08 ± 0.04 ± 0.02 


0.95 


57+ u - ia 

U - 0/ -0.14 



1 Introduction 

The hadronic B decays have been studied for many years since they offer an excellent place to study 
the CP violation and search for new physics hints [1]. The hadronization of the final states is non- 
perturbative in nature, and the essential problem in handling the decay processes is the separation of 
different energy scales, namely, the so-called factorization assumption. Many factorization approaches 
have been developed to calculate the B meson decays, such as the naive factorization [2], the general- 
ized factorization [3, 4], the QCD factorization [5], as well as the perturbative QCD approach (PQCD) 
based on kx factorization [6, 7]. Most factorization approaches are based on heavy quark expansion 
and light-cone expansion, only the leading power or part of next to leading power contributions are 
calculated to compare with the experiments. Nevertheless for the penguin-dominated decay channels, 
the power corrections and the nonperturbative contributions may be large, since the theoretical pre- 
dictions for some channels cannot fit the data quite well. There are some problems, such as the sin 2(3 
problem in penguin dominated modes [8], which suggests that more dynamics of penguin dominating 
B decays should be studied. 

Recently, with more and more data, the B factories have measured some decays that the final 
state contains two vector mesons [9, 10, 11]. In the B — > VV modes, both the longitudinal and 
the transverse polarization can contribute to the decay width, and the polarization fractions can be 
measured by the experiments. The naive counting rules based on the factorization approaches predict 
that the longitudinal polarization dominates the decay ratios and the transverse polarizations are 
suppressed [12] due to the helicity flips of the quark in the final state hadrons. But some data shown 
in table 1 are quite different from the theoretical predictions for the penguin dominated modes. 

The small longitudinal polarization fraction in B — > 4>K* decays has been considered as a puzzle, 
many theoretical efforts have been performed to explain it [13, 14, 15, 16, 17, 18, 19, 20]. In PQCD 
approach, the coefficients of penguin operators have been evolved to the scale of about v/AMg, so 
these coefficients become larger compared to the factorization approach, in which the hard scale are 
at the scale of Mb, so that the penguins' contribution are enhanced in PQCD approach. Besides, 
the annihilation diagrams, which is power suppressed in QCD factorization, are also included. Thus 
the PQCD approach can give a larger branching ratio and fits the experiments well in B — ► PP, PV 
case. For B — ► (f)K*, the annihilation diagram with the (S + P)(S — P) type operators will break the 
naive counting rules [15], the transverse polarization is enhanced to about 0.25. But the branching 
ratios calculated in the PQCD approach [21] are too large if we adopt the old K* meson's parton 
distribution amplitudes derived from QCD sum rules. As mentioned in [22], things will get better (59% 
of longitudinal polarization) if we adopt the asymptotic form of the K* meson's parton distribution 
amplitudes. 

In this paper, we will perform the leading order PQCD calculation of penguin dominated processes 
B — > pK* and B — ► ujK* . The branching ratios have been measured by the B factories [23] which are 
given in table 2. And the measured CP asymmetries are: Acp{B + — > p + K*°) = —0.14 ± 0.17 ± 0.4 
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Table 2: Branching ratios (10 6 ) of B — > * measured by the B factories 



Process 


BaBar 


Belle 


world average 


5 U -► 


< 24 




< 24 


5+ p u K*+ 


10.6ti^ ± 2.4 




10.6lij;jj 


5+ -► p+K*° 


17.0 ± 2.9i^ 


8.9 ± 1.7 ± 1.0 


10.5 ± 1.8 






< 2.6 


< 2.6 


5+ -► ujK*+ 


< 7.4 




< 7.4 


b° -► wif* u 


< 6.0 




< 6.0 



and A C p(B+ p°K* + ) = 0.20t° H ± 0.04. These channels have been studied within the QCD 
factorization framework, but the predictions are not quite consistent with the data, especially the 
polarization fractions [24]. We hope the PQCD approach could give a better theoretical prediction. 

The paper is organized as follows: In Sect. 2 we will present the framework for three scale PQCD 
factorization theorem. Next we will give the perturbative calculation result for the hard part. In Sect. 
4, numerical calculation for branching ratio and CP violation are given. Final section is devoted to 
summary. 



2 The Theoretical Framework 

The PQCD factorization theorem has been developed for non-leptonic heavy meson decays [25], based 
on the formalism by Brodsky and Lepage [26], and Botts and Sterman [27]. In the two body hadronic 
B decays, the B meson is heavy, sitting at rest. It decays into two light mesons with large momenta. 
Therefore the light mesons are moving very fast in the rest frame of B meson. 

To form the fast moving final state light meson, in which the two valence quarks should be collinear, 
there must be a hard gluon to kick off the light spectator quark d or u in the B meson (at rest). So the 
contribution from the hard gluon exchange between the spectator quark and the quarks which form 
the four quark operator dominates the matrix element of the four quark operator between hadron 
states. This process can be calculated perturbatively, but the endpoint singularity will appear if 
we drop the transverse momentum carried by the quarks. After introducing the parton's transverse 
momentum, the singularity is regularized, and additional energy scale is present in the theory, then the 
perturbative calculation will produce large double logarithm terms, these terms are then resummed 
to the Sudakov form factor. The uncancelled soft and collinear divergence should be absorbed into 
the definition of the meson's wave functions, then the decay amplitude is infrared safe and can be 
factorized as the following formalism: 



C(t) x H(t) x x exp 



-s(P,b)-2 f % q (a s (fl)) 



(1) 



where C(t) are the corresponding Wilson coefficients of four quark operators, are the meson 

wave functions and the variable t denotes the largest energy scale of hard process H, it is the typical 
energy scale in PQCD approach and the Wilson coefficients are evolved to this scale. The exponential 
of S function is the so-called Sudakov form factors, which can suppress the contribution from the 
nonperturbative region, making the perturbative region give the dominated contribution. The "x" 
here denotes convolution, i.e., the integral on the momentum fractions and the transverse intervals of 
the corresponding mesons. Since logarithm corrections have been summed by renormalization group 
equations, the factorization above formula does not depend on the renormalization scale p explicitly. 
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In the resummation procedures, the B meson is treated as a heavy-light system. In general, the 
B meson light-cone matrix element can be decomposed as [28, 5] 



L 



1 —^(Q\b a {<d)dp(z)\B(j) B )) 



(2vr) 4 

[{£>B + m B )ls [fo(ki)- A+ ^- ^ B (k!)]} , (2) 



where n+ = (1, 0, Ot), and n_ = (0,1, Ot) are the unit vectors pointing to the plus and minus direc- 
tions, respectively. As pointed out in ref.[29], this kind of definition will provide light-cone divergence, 
and more involved studies have been performed [30, 31]. Here we only use it phenomenologically to 
fit the data, so we still use the old form. From the above equation, one can see that there are two 
Lorentz structures in the B meson distribution amplitudes. They obey the following normalization 
conditions 

In general, one should consider both of these two Lorentz structures in calculations of B meson decays. 
However, it can be argued that the contribution of <j>B is numerically small [32, 33], thus its contribution 
can be neglected. Therefore, we only consider the contribution of Lorentz structure 

in our calculation. Note that we use the same distribution function (f>B(ki) for the fas term and 
the tub term in the heavy quark limit. For the hard part calculations in the next section, we use 
the approximation rrit, ~ tub, which is the same order approximation neglecting higher twist of 
(tub — m>b) /tub- Throughout this paper, we take light-cone coordinates, then the four momentum 
P^ = (p° = ' = P 3 )/v / 2) and pt = (p 1 ,p 2 )- We consider the B meson at rest, the momentum is pb = 
(ms/\/2)(l, 1,0t)- The momentum of the light valence quark is written as (kf , fcj~,kir), where the 
kix is a small transverse momentum. It is difficult to define the function 4>B(kf, k± , kir). However, 
the hard part isn't always dependent on kf if we make some approximations. This means that kf 
can be simply integrated out for the function 4>s(kf , \<.it) as 

(f>B(x,k 1T ) = J ' dkf(j) B {kf ,k± ,ki T ) (5) 

where x = /p^ is the momentum fraction. Therefore, in the perturbative calculations, we do not 
need the information of all four momentum k±. The integration above can be done only when the 
hard part of the subprocess is independent on the variable kf. 

The K* and p mesons are treated as a light-light system. At the B meson rest frame, they are 
moving very fast. We define the momentum of the K* as P2 = (ms/V2)(l — r|,r|,0y). The p 
has momentum P 3 = (tub /V2) (r| , l — r^Or), with r 2 = Mk*/Mb and r 3 = M^/Mb- The light 
spectator quark in K* meson has a momentum (kf , 0, \s~2t)- The momentum of the other valence 
quark in this final meson is thus (P 2 + — kf,0, — k 2 y). The longitudinal polarization vectors of the K* 
and p are given as: 

P 2 Mk* m P3 M p 

e2{L) = — — n_ , e 3 (L) = — - — n+ , (6) 

M K * P 2 ■ n_ M p P 3 • n + 

which satisfy the normalization e\(L) = e| (L) = — 1 and the orthogonality e 2 (L) ■ P 2 = e^(L) ■ P3 = 
for the on-shell conditions P 2 = and Pf = Mi 2 . We first keep the full dependence on the light 
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meson masses Mk* and M p with the momenta Pi and P3. After deriving the factorization formulas, 
which are well-defined in the limit Mx*,M p — > 0, we drop the terms proportional to r 2 ,, r 2 ^* ~ 0.04. 
The transverse polarization vectors can be adapted directly as 

= ^(0, 0, 1, i) , e(-) = i=(0, 0, 1, -i) . (7) 

If the K* meson (so as to other vector mesons) is longitudinally polarized, we can write its wave 
function in longitudinal polarization [32, 34] 

<K*-(P,e L )\s a (z)up(0)\0> 

= ~^=^ J q dxe lxP ' z y^K-^K-ix) Imjf.^.(i)] + m K *<t>K*( x )} ■ ( 8 ) 

The second term in the above equation is the leading twist wave function (twist-2), while the first and 
third terms are sub-leading twist (twist-3) wave functions. If the K* meson is transversely polarized, 
its wave function is then 

<K*-(P,e T )\s a (z)up(0)\0> = dxe txP - z {^[^<p T K ,(x)+m K ^ v K ,(x) 

+im K *e lxvpa w fX e u n p v a (f) a K * (x)} . (9) 

Here the leading twist wave function for the transversely polarized K* meson is the first term which 
is proportional to . 

The transverse momentum kjT is usually converted to the b parameter by Fourier transformation. 
The initial conditions of <fii(x), i = B, K*, p, are of nonperturbative origin, satisfying the normalization 

\ i (x,b = 0)dx = ^=f i , (10) 

with fi the meson decay constants. 



3 Perturbative Calculations 

With the above brief discussion, the only thing left is to compute the hard part H . We use the 
notation M\ =< Vi{\)V2{X)\H^ \B > for the helicity matrix element, A = 0, ±1. For decays of B to 
two vector mesons, the amplitude can be expressed by three invariant helicity amplitudes, defined by 
the decomposition 

M A = M«e^(A) • e ;(A) + M^e* K «{\) • P p e* p (\) • P K * + M^tE pv ^{\)e™ {\)P£«P° . (11) 

According to the naive counting rules mentioned before, we can estimate that polarization fractions 
satisfy the relation: |Mo| 2 3> |M_| 2 3> |M + | 2 . These three helicity amplitudes can be expressed as 
another set of helicity amplitudes, 

M = M 2 B M L , Af± = M^MntM^V? - 1M T (12) 

where the Ml, Mn and Mp can be extracted directly from calculation of the Feynman diagrams, and 
r' = ■ The formula for the decay width is 

r = dbs2X>"w- < 13 > 
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Here p is the absolute value of the 3-momentum of the final state mesons. And we have 

E M (t) M M = |M | 2 + |M + | 2 + |M_| 2 .. 
The weak Hamiltonian H e g for the AB = 1 transitions at the scale smaller than m\v is given as 

[35] 

/ 10 \ 
VvbKa (CiOi + C 2 2 ") - V tb v: s [J2CiOi + c g o g 



G F 



(14) 



We specify below the operators in H e jf for b — > s: 



Of = Sa^Lup ■ , 0% = s a ^Lu a ■ u^^Lbp , 

03 = Sa^Lba ■ Eg' Qpl/J-Lq'fS , O4 = SaJ^Lbp ■ Eg' Qpl/xWa > 

O5 = s a ^Lb a ■ Eg' q'piixRq'p , Oq = sa^Lbp ■ Eg' q'pinRq'a > (15) 

7 = ^Sa^Lba • Eg' Q'plnWp ' °8 = \s a ^Lbfi ■ Eg' e^q'^^Rq^ , 
O9 = §s Q 7^a • Eg' eq'q'pl^Lqp , O 10 = paj^Lbp ■ Eg' e^q'^^Lq^ . 

Here a and /3 are the SU(3) color indices; L and i? are the left- and right-handed projection operators 
with L = (1 — 75), R = (I + 75). The sum over q' runs over the quark fields that are active at the 
scale n = 0(mb), i.e., (q'e{u, d, s, c, b}). 

The diagrams for these decays are completely the same as ones in the decay B — > Kir. Here we 
take the decay B — > p°K* + as an example, whose diagrams are shown in Figure 1. These are all 
single hard gluon exchange diagrams, containing all leading order PQCD contributions. The analytic 
calculation is performed through the contraction of these hard diagrams and the Lorenz structures of 
the mesons' wave functions. The first row and the third row in Figure 1 are called emission diagrams, 
with the p meson or K* meson emitted. The analytic formulae for the K* meson emission diagram 
is exactly the same as the emission diagrams of B — ► K* + (f> with — > fx*, fx* — > fp, and we can 
get the formulae for the p emission diagrams through the change — ► f p , X3 — > x 2 from B — ► K* + cp. 
As to the annihilation diagrams, we make the same change as the if* emission diagrams for the 
corresponding diagrams of B — ► i"C* + <^>, then we can get the right analytic formulae. 

In PQCD approach, only Wilson coefficients are channel dependent. There are six different decay 
channels in B + (B°) — > p(oj)K* decays, and the B~(B°) decays are their CP conjugation. All these 
decays are included in the twelve diagrams, the only changes needed are external quarks and the 
Wilson coefficients. We summarize the Wilson coefficients for each channels in table 3. In this table 
the coefficients are defined as 

ai = Ci + C 2 /N c , a 2 = C 2 + d/N c , 
4 = C q 3 + C q JN c + ^e q (C 9 + C w /N c ), 

a\ = Cl + Ct/N c + ^e q (C w + C 9 /N c ), 

4 = Cl + Cl/N c + ^e q (C 7 + C 8 /N c ), 

4 = Cl + C q 5 /N c + ^e q (C 8 + C 7 /N c ), (16) 



and 



a 3 q = C$ + ^e q C 9 , (17) 
3 

a'i = Cf + -e q C 10 , 
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(a) (b) (c) (d) 




(f) (g) (h) 




(i) (j) (k) (1) 



Figure 1: Diagrams contributing to the B + — » K* + p decays 



Table 3: Wilson coefficients (the characters in the first row stand for the diagrams in figure 1 ) 



Process 


(a)(b) 


(c)(d) 


(g)(h) 


(e)(f) 


(00) 


000) 


B° -» p-K*+ 






(d) (d) 


af ,a? J 


(«) 

a 2 , a 4 




B + -> p°K*+ 


ai,a 3 , a 5 


c 2 ,o2 u - d) ,o; (u - dJ 




Oi,a 3 ,a 5 


a 2 , a 4 


Oi,a 3 ,a 5 


B + -> 






(it) (u) 
a2,a 4 ,a 6 


ci,ar,ar 


a? 






01,0^ , Og 


r< A u ~ d ) A u - d ) 


o (d) a W 








5° -► ujK*° 


„ „(«+<*) 

Ol , Og , Og 


O2,o 4 ,a 6 


(d) (d) 
a 4 > a 6 


a V) a' (d} 




'(<*) '(^) 


B + -> 


ai,a 3 ,a 5 


o 2 ,a 4 ,a 6 


n >> „W 

a2,a 4 ,a 6 


Ci,a 3 ,a 5 


(u) 
«2, 04 


Oi,a 3 ,a 5 



ag = CZ + ^e q C 8 . (18) 
4 Numerical Calculations and Discussions of Results 

In the numerical calculations we use [36] 

f B = WOMeV, m K * = 0.892 GeV, m p = 0.77 GeV, 
M B = 5.28GW, f K » = 217MeV, /£. = WOMeV, 
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M w = 80AlGeV, 
m w = 0.782 GeV, 



f p = 205MeV, 
f u = 195MeV, 



t b ± = 1.671 x 10 



-12 



T B 



1.536 x 10- 12 s, A& 4) 



ft = !55MeV, 
fl = UOMeV, 

250MeV. 



MS 



(19) 



The distribution amplitudes (p p (x) (0^(x)) and (p l K *(x) of the light mesons used in the numerical 
calculation are listed in Appendix A. 

For B meson, the wave function is chosen as 



4>B(x,b) 



N B x 2 (l - x) 2 exp 



Ml x 2 
2u,, 2 



1 



(u b bf 



(20) 



with uj b = 0.4 GeV [37], and the normalization constant N B = 91.784 GeV. We would like to point 
out that the choice of the meson wave functions and the parameters above is the result of a global 
fitting for B — ► irir and B — > ttK decays [6, 7]. 

For the CKM matrix elements, we use \V us V* b \ = 0.00078, \V ts V t l\ = 0.0395. We leave the CKM 
angle 03 as a free parameter, which is defined as 



V ub = \V ub \exp(-i(p 3 ). 
The decay amplitude of B — > K*p can be written as 

-M (i) = V,: h V us T^ -V t * b V ts P^ 



(21) 



'ub v us- 

V*,V 

v ub v us ± 



1 



z (i) e i(-<t>3+S (l) ) 



i = 1,2,3, 



(22) 



where = y^^y" y$) ' an< ^ ^ ls ^ e re l & ti ve strong phase between tree (T) diagrams and penguin 
diagrams (P). and <jW can be calculated perturbatively. Here in PQCD approach, the strong phases 
come from the non-factorizable diagrams and annihilation type diagrams (see (c) ~ (h) in Figure 1). 
The internal quarks and gluons can be on mass shell, and then poles appear in the propagators, which 
can provide the strong phases. The predominant contribution to the relative strong phase 5 comes 
from the annihilation diagrams, (g) and (h) in Figure 1. 

This mechanism of producing strong phase is very different from the so-called Bander-Silverman- 
Soni (BSS) mechanism [38], where the strong phase comes from the perturbative charm penguin 
diagrams. The contribution of BSS mechanism to the direct CP violation in B — > K*p is only in the 
higher order corrections (a s suppressed) in our PQCD approach. Therefore we can safely neglect this 
contribution. 

The corresponding charge conjugate B decay is 

M W = V ub V: s T^ - V tb V t * s PW 
= V ub V*T^\l 



(23) 



In contrast to the decay of B to pseudoscalar mesons like B — > Kir, where the decay widths can be 
expressed in terms of 5 and 03 in a simple way, here for B decay to two vector mesons, there are 3 
types of amplitudes, and this makes the dependence of decay widths on 5 and 03 very complicated. 
The averaged decay width for B and its CP conjugation decays can be expressed as a function of a 
CKM phase angle 03. 



r = V^72 \ V ub V us\ 2 [Tl(l + zl + 2z L cos<j )3 cos 5 L ) + 2 V T?{1 + z\ + 2z iCOS 3 cos 5 t )\. 

87TM B i^T 



(24) 
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Figure 2: Averaged branching ratios (10 6 ) of B — > K* p and i? — > as a function of CKM angle 
03, where the lines B,C,E,F in diagram (a) represent B + -> p+K*°,B + -> B° -> p ^* , 

5° — > p~K* + respectively, and in diagram (b), B denotes 5° — > wiC* , C denotes -B + — > 

From this formula we can know that when contribution from the penguin diagrams are much larger 
than that from the tree diagrams, i.e., S> l,i = L,N,T, then the branching ratios are insensitive 
to the angle, but when they are comparable, the dependence on 03 will be strong. We show the 
branching ratios of these decays in Figure 2, from which we can see that the penguin dominant decays 
B + — ► p + K*°,B° — ► p°K*°, and — > u;.ftr* are almost independent on 03, but the dependence on 
03 of the other three channels is strong, because of the tree and penguin interference. 

In Figure 3, we plot the dependence of longitudinal polarization fractions Tl/T on the CKM angle 
03. We find that this quantity is not very sensitive to 03 in all decay channels. If we fix 03 at about 
60°, we find that for the decays B + — > K* + p° and B + — > K*°p + , the longitudinal fractions are 0.89 
and 0.82 respectively. As mentioned before, we calculate the annihilation type diagrams in PQCD 
approach. If the four quark operator has the Dirac structure like (S — P)(S + P), there is no helicity 
flip suppression to the transverse polarization, so that the longitudinal fractions are considerably 
suppressed. One can see that our results for B + — > K*°p + are consistent with BaBar, but different 
from Belle (we hope more efforts from experimental side to test our prediction). As to B + — ► K* + p°, 
our result is a little smaller, but still agree with the data within the la error bar. 

The new analysis of the K* meson wave function from QCD sum rules [39] shows that the leading 
twist distribution amplitude 4>k* (x) of longitudinal polarization should be very close to the asymptotic 
one. According to Li's suggestion [22], we test our result using the asymptotic wave functions for the 
longitudinal polarization part. The numerical results are given in table 4. We find that the longitudinal 
fraction and the branching ratios for all the channels are reduced. Note that Figure 2 shows that the 
branching ratios of B° — > K*°uj and B + — > K* + oj are larger than the experimental limits where we 
use the wave functions given in the appendix. But if we adopt the asymptotic form, the branching 
ratios decrease. Comparing the table with the experimental data, it seems that the asymptotic form 
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Figure 3: Longitudinal polarization fraction of B — > K*p and B — ► iT*u; as a function of CKM angle 

03 



Table 4: Branching ratios (10 6 ) and polarization fractions using different type of light meson wave 
functions (the CKM phase angle 4>3 is fixed as 60°; the w.f. stands for wave function) 



Quantity 


w.f. in the appendix 


asymptotic w.f. 


Br(B u -► p~K* + ) 


13 


9.8 


Br(B+ -► p+ivr* u ) 


17 


13 


Br(B+ -► p u £T*+) 


9.0 


6.4 


5r(B u -► p u K* u ) 


5.9 


4.7 


5r(B+ -► ^*+) 


7.9 


5.5 


5r(B u -► c^* u ) 


9.6 


6.6 


i? L (5 u p-ir<+) 


0.78 


0.71 


R L (B+ -► p+if* u ) 


0.82 


0.76 


i? L (S+ -► p u if*+) 


0.85 


0.78 


i? L (5+ -► 


0.81 


0.73 


i? L (5 u p u K* u ) 


0.74 


0.68 


i?i(5 u wiT* u ) 


0.82 


0.74 



is more convincing. More study of the vector meson's wave functions are required. 

It has been confirmed that there is big direct CP violation in B — ► ttK and B — > irir decays [23], 
and the PQCD approach can give right predictions from the annihilation topology [40] rather than 
the BSS mechanism. Here we take the definition (note that our definition has opposite sign when 
comparing with the definition used in [23]) 



A C p = 



T(B -> f) - T(B -> f) 
F(B — > /) + T(B — »■ /) 



(25) 
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The direct CP violation parameters as a function of (fo are shown in Figure 4. Since CP asymmetry is 
sensitive to many parameters, the line should be broadened by uncertainties. The direct CP violation 
parameter of B + — > K* + p°, B° — > K* + p~ and 2? + — > i^T* + cj can be large as 15 - 20% when <p 3 is 
near 60°, but for B + — > K*°p + , the direct CP violation is very small for the very tiny tree diagram 
contribution. The final state is not the CP eigenstate, so the mixing induced CP violation is more 
complicated, and we do not give it here. 



0.10 r 

0.08 ■ 

0.06 ■ 

0.04 ■ 

0.02 • 
0.00 

-0.02 ■ 

-0.04 ■ 

-0.06 ■ 

-0.08 • 
-0.10 
-0.12 

-0.14 • 

-0.16 ■ 
-0.18 

-0.20 ■ 



-pK" 
pK* 




-ojK* 
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Figure 4: Direct CP violation of B — > K*p and B — > K*uj as a function of CKM angle 

The angular distributions depend on the spins of the decay products of the decay vector mesons 
K* and p. For example, for B + — ► K* + p° — > (i^7r)(7r + 7r _ ) the differential decay distribution is [41] 

d 3 r 9 rir T . 2 . 2 r L 2 2 

— 1 < -— ■ sm U\ sm t/2 + -pr " cos #1 cos 6/2 



dcosMcos0 2 # 8vr 14 T ' T 

+ - sin sin 26*2 [oti ■ cos <f> — fii ■ sin 0] 

+i sin 2 6»i sin 2 6 2 [a 2 ■ cos 20 - ■ sin 2</> ] J . (26) 

In (26) 6\ is the polar angle of the K in the rest system of the K* with respect to the helicity axis. 
Similarly 2 is the polar angle of the tt + in the p rest system with respect to the helicity axis of 
the p°, and (j> is the angle between the planes of the two decays K*~ — ► Kit and p° — > 7r + 7r~. The 
coefficients in the decay distribution are related to the helicity matrix elements by 

rr [M+xp + IM-xl 2 _ |M | 2 

T |Mo| 2 + |M+i| 2 + |M_iP' T \Mo[ I + \M +1 \^ + \M- 1 [^ 

Re(M +1 M*+M. 1 M*) Im{M +l M*-M^M*) , . 

ttl ~~ IMo^ + IM+iP + IM-iP' ^ _ |M p + |M+iP + |M_iP ' y Z ') 

Re(M +1 M*_ 1 ) R ImjM+iM^) 

012 ~ |M |^ + |M + i| 2 + |M_iP' ^ 2 - |M |^ + |M +1 p + |M_i|^ 
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The integration over angles 9\, 62 in (26) yields the (f> distribution of the decay width 

2vr— = T(l + 2a 2 cos2(/>- 2/3 2 sin20), (28) 

where the coefficients a 2 , 02 can be obtained from (27) by using the M\ which is calculated in PQCD 
approach. Because a 2 and 02 are very small, the decay width is almost independent in (ft, then the 
CP violation from the angular distribution will be very tiny in the standard model. 



5 Summary 

We performed the calculations of B + -► K*+p°, B + -► K*°p+, B° -► K*+p~, B° -► K*°p° and 
B + — > ivT* + oj, f?° — ► ivT*°aj in PQCD approach. In this approach, we calculated the non-factorizable 
contributions and annihilation type contributions in addition to the usual factorizable contributions. 

We found that the annihilation contributions were not so small as expected in a simple argument. 
The annihilation diagram, which provides the dominant strong phases, plays an important role in 
the direct CP violations. We expect large direct CP asymmetry in the decays of B + — > K* + lj, 
B° — > K* + p~ and B + — ► K* + p°. We also study the helicity structure and angular distribution of the 
decay products. The current running B factories in KEK and SLAC will be able to test the theory. 
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A Wave Functions of Light Mesons Used in the Numerical Calcu- 
lation 

For the light meson wave function, we neglect the b dependence part, which is not important in 
numerical analysis. We choose the different distribution amplitudes of p meson longitudinal wave 
function as [34], 



4> P (x) = 6/ p x(l - 

= fJ{st 2 + o. 



l + 0.18C 2 3/2 (t) 



3t 2 



5t 2 -3 



+ 0.21 



3 - 30t 2 + 35t 4 ] } , 



(j) s p ( x ) = 3/J t [l + 0.76(10x 2 - 10s + 1)] , 
where t = 1 — 2x. The Gegenbauer polynomials are defined by 



d /2 (t) = i(3t 2 - 1), C\' 2 {t) = |(35t 4 - 30t 2 + 3), 



Cl ,2 (i) 



|(5t 2 - 1), Cl ,2 {t) = f(21i 4 - 14t 2 + 1). 



For the transverse p meson we use [34]: 



= 6/Jx(l-x) 



1 + 0.2C 2 3/2 (t) 



f p |^(1 + t 2 ) + 0.24(3t 2 - 1) + 0.12(3 - 30t 2 + 35t 4 )| 
1 + 0.93(10x 2 - lOx + 1) 



(29) 
(30) 
(31) 

(32) 

(33) 
(34) 

(35) 
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For the u meson, we use the same as the above p meson, except changing the decay constant f p with 

foj- 

We choose the light cone distribution amplitudes of K* meson longitudinal wave function as [34], 

4> K *{x) = 6f K *x(l -x) [l + 0.57t + 0.07C2 /2 (t)l , (36) 

tf K «(x) = /£*{o.3/(3t 2 + 10/-l) + 1.68C 4 1/2 (t) + 0.06t 2 (5t 2 -3) 

+0.36[l-2t-2tln(l-x)]}, (37) 

cf) s K , ( x ) = f£, {3/ [l + 0.2/ + 0.6(Wx 2 - lQx + 1)] - 0.12x(l - x) 

+0.36[1 - 6x - 2 ln(l - x)]} . (38) 

The light cone distribution amplitudes of K* transverse wave function are used as [34] 

4> T K *{x) = 6/£.x(l -x) [l + 0.6/ + 0.04(72 /2 (t)l , (39) 

<t>K*{x) = /^|^(l + / 2 + 0.44/ 3 )+0.4C 2 1/2 (/)+0.88C 4 1/2 (/) 

+0.48[2x + ln(l -x)]}, (40) 

<f> a K .{x) = ^{3t[l + 0.19/ + 0.81(10x 2 -10x + l)j -l.Ux{l-x) 

+0.48[1 - 6x - 2 ln(l - x)]} . (41) 
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